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Abstract 

We introduce a vertex amplitude for 4d loop quantum gravity. We derive it from a conventional 
quantization of a Regge discretization of euclidean general relativity. This yields a spinfoam sum 
that corrects some difficulties of the Barrett-Crane theory. The second class simplicity constraints 
are imposed weakly, and not strongly as in Barrett-Crane theory. Thanks to a flip in the quantum 
algebra, the boundary states turn out to match those of SO(3) loop quantum gravity - the two can 
be identified as eigenstates of the same physical quantities - providing a solution to the problem of 
connecting the covariant SO(<i) spinfoam formalism with the canonical SO(3) spin-network one. 
The vertex amplitude is SO(3) and S l O(4)-covariant. It rectifies the triviality of the intertwiner 
dependence of the Barrett-Crane vertex, which is responsible for its failure to yield the correct 
propagator tensorial structure. The construction provides also an independent derivation of the 
kinematics of loop quantum gravity and of the result that geometry is quantized. 

1 Introduction 

While the kinematics of loop quantum gravity (LQG) [T] is rather well understood [2J [3] , its dynamics 
is not understood as cleanly. Dynamics is studied along two lines: hamiltonian [3] or covariant. The 
key object that defines the dynamics in the covariant language is the vertex amplitude, like the vertex 
amplitude = ej^S(pi + P2 + k) that defines the dynamics of perturbative QED. What is the 

vertex of LQG? 

The spinfoam formalism [31 O [71 [HI can be viewed as a tool for answering this question: the 
spinfoam vertex plays a role similar to the vertices of Feynman's covariant quantum field theory. 
This picture is nicely implemented in three dimensions (3d) by the Ponzano- Regge model |10j , whose 
boundary states match those of LQG [11] and whose vertex amplitude can be obtained as a matrix 
element of the hamiltonian of 3d LQG [T^]. But the picture has never been fully implemented in 
4d. The best studied model in the 4d euclidean context is the Barrett-Crane (BC) theory [8], which 
is based on the vertex amplitude introduced by Barrett and Crane [7J. This is simple and elegant, 
has remarkable finiteness properties [13] . but the suspicion that something is wrong with it has long 
been agitated. Its boundary state space is similar to, but does not exactly match, that of LQG; in 
particular the volume operator is ill-defined. Worse, recent results [2] indicate that it appear to fail 
to yield the correct tensorial structure of the graviton propagator in the low-energy limit [T3] • 

It is then natural to try to correct the BC model [T51 [TTJ [T5] . The difficulties are all related to the 
fact that in the BC model the intertwiner quantum numbers are fully constrained. This follows from 
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the fact that the simplicity constraints are imposed as strong operator equations (C n ip = 0). However, 
these constraints are second class and it is well known that imposing second class constraints strongly 
may lead to the incorrect elimination of physical degrees of freedom [19 . In this paper we show that 
the simplicity constraints can be imposed weakly ((<fi C' n ip) = 0), and that the resulting theory has 
remarkable features. First, its boundary quantum state space matches exactly the one of 5*0(3) LQG: 
no degrees of freedom are lost. Second, as the degrees of freedom missing in BC are recovered, the 
vertex may yield the correct low-energy n-point functions. Third, the vertex can be seen as a vertex 
over 50(3) spin networks or 50(4) spin networks, and is both 5*0(3) and SO (A) covariant. 

These results have been anticipated in a letter [20] . Here we derive them via a proper quantization 
of a discretization of euclidean general relativity (GR). Indeed, although spinfoam models have been 
derived in a number of different manners [H [3 12 lj , most derivations involve peculiar procedures 
or intuitive and ad hoc steps. It is hard to find a proper derivation of a spinfoam model from the 
classical field theory, which follows well-tested quantization procedures. Here we try to fill this gap. 

From the experience with QCD, one can derive the persuasion that a nontrivial quantum field 
theory should be related to a natural lattice discretization of the corresponding classical field theory 
[2"2"] . This persuasion is reinforced by the LQG prediction of an actual physical discretization of 
spacetime. Here, we reconstruct the basis of the euclidean spinfoam formalism as a proper quantization 
of a lattice discretization of GR. Conventional lattice formalisms as the ones used in QCD are not very 
natural for GR, since they presuppose a background metric. Regge has found a particularly natural 
way to discretize GR on a lattice [23], known as Regge calculus. Quantization of Regge calculus has 
been considered in the past [24] and its relation to spinfoam theory has been pointed out (see [25] 
and references therein). Here we express Regge calculus in terms of the elementary fields used in 
the loop and spinfoam approach, namely holonomies and the Plebanski two-form, and we study the 
quantization of the resulting discrete theory (on lattice derivations of loop gravity, see [6l [27l [28] ) . 

The technical ingredient that allows a nontrivial intertwiner state space to emerge can be inter- 
preted as a "flip" of the 50(4) algebra, namely an opposite choice of sign in one of its two SU (2) 
factors. The possibility and the relevance of a "flipped" symplectic structure was noticed by Baez 
and Barrett in [26 (where they attribute the observation to Jose- Antonio Zapata) and by Montesinos 
[29] , Montesinos, in particular, has discussed the classical indeterminacy of the symplectic structure 
in detail. The flip can be viewed as the equivalent of the Ashtekar "trick" , which yields a connection 
as phase space variable. The 50(4) generators turn out to directly correspond to the bivectors asso- 
ciated to the Regge triangles, rather than to their dual. Using this, we find a nontrivial subspace of 
the SO (4) intertwiner space, which corresponds to closed tetrahedra and maps naturally to an 50(3) 
intertwiner space. 

This path leads to a quantum theory that appears to improve several aspects of the better known 
spinfoam models. In particular: (i) the geometrical interpretation for all the variables becomes fully 
transparent; (ii) the boundary states fully capture the gravitational field boundary variables; and 
(iii) correspond precisely to the spin network states of LQG. The identification is not arbitrary: the 
boundary states of the model are precisely eigenstates of the same quantities as the corresponding 
LQG states. This last result provides a solution to the long-standing difficulty of connecting the 
covariant 50(4) spinfoam formalism with the 50(3) canonical LQG one. It also provides a novel 
independent derivation of the LQG kinematics, and, in particular, of the quantization of area and 
volume. Finally, (iv) the vertex of the theory is similar to, but different from, the BC vertex, leading 
to a dynamics that might be better behaved in the low-energy limit. 

The paper is organized pedagogically and is largely self contained. Section 2 reviews background 
material: properties of 50(4) and its selfdual/anti-selfdual split, and the definition of the fields and 
the formulation of classical GR as a constrained Plebanski theory. In section 3 we discretize the 
theory on a fixed triangulation of spacetime. We do so simply by taking standard Regge calculus and 
re-expressing it in terms of the (discretized) Plebanski two-form. The resulting theory is governed 
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by the geometry of a 4-simplex, which we illustrate in detail. All basic relations among the variables 
have a simple interpretation in these terms. The 10 components of the metric tensor g a b in a point 
(or in a cell) can be interpreted as a way to code the 10 variables determining the shape of the cell. 
In particular, the norms of the discretized B fields on the faces are their areas and the scalar product 
on adjacent triangles codes the angle between the triangles. While these "angles" and "areas" are 
independent in BF theory, they are related if they derive from a common metric, namely in GR. 
In Section 4 we study the quantization of the system. We explain the difficulties of imposing the 
constraints strongly, study the weak constraints and write their solution. Finally we construct the 
vertex amplitude. 

We work in the euclidean signature, and on a fixed triangulation. The issues raised by recovering 
triangulation independence and the relation with the Lorentzian-signature theory will be discussed 
elsewhere. 



2 Preliminaries: Plebanski two-form and structure of 5*0(4) 

Riemannian general relativity (GR) is defined by a riemannian metric g a b(x), where a,b — 1,2,3,4 
and the Einstein-Hilbert action 

S[g] = J^R = J ^g ab R ab . (1) 

where g ab is the inverse, g the determinant, and R a b the Ricci curvature of g a t,. A good number of 
reasons, such as for instance the fact that this metric formulation is incompatible with the coupling 
with fermions, suggest to use the tetrad field ej(x), I = 1, 2, 3, 0, (the value instead of 4 is for later 
convenience and does not indicate a Lorentzian metric) or its inverse, namely the tetrad one-form field 
e 1 (x) — e I a (x)dx a , to describe the gravitational field. This is related to the metric by e^e^ = g a b- Sum 
over repeated indices is understood, and the up or down position of the I indices is irrelevant. The 
spin connection of the tetrad field is an SO (4) connection w / ' 7 [e] satisfying the torsion- free condition 

De 1 = de 1 + lu 1 j[e] A e J = 0. (2) 

The action |T]) can be rewritten as a function of e 1 in the form 

S[e] = J (det e) e?e^F#[w[e]] = \ J e IJKL e 1 A e J A F KL [u[e]] (3) 

where F[uj] is the curvature of to. Alternatively, GR can be defined in first order form in terms of 
independent variables uj ij and e 1 , by the action 

S[e,w] = \Je U KL e 1 Ae J AF KL [u}. (4) 

In this case, ([2]) is obtained as the equation of motion for u. 



2.1 Plebanski two-form and simplicity constraints 

2^ab 



At the basis of the spinfoam formalism is the use of the Plebanski two-form E /J = h^ii dx a A dx b , 



defined as 

S /J — e 1 A e J . (5) 

or its dual, usually called B IJ = \B 1 J ) dx a A dx b for a reason that will be clear in a moment, defined 

as 

B = -e kl — ^ e e A e yo) 
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We use the following notation for two-index objects: a scalar product: A ■ B = A IJ Bjj; a norm: 
\B\ 2 = B ■ B, and the duality operation (*B) IJ = \e u ' KL B KL . So that, in particular, *B ■ B = 
\tUKL B IJ B KL . Thus we write © in the form 

B = *E = *(e 7 Ae J ) (7) 



The geometrical interpretation of the Plebanski two-form (or the B two-form) is captured by the 
following. Observe that 

I S afc 1 2 = \B ab \ 2 = g aa 9bb - QabQab = 2A 2 fc (8) 

and 

E a h ' ^ac = B a b ' B ac — QaaQbc ~ 9ab9ac = 2Jaabc (9) 

The quantity A a \, gives the area element A a \ > dx a dy h of the infinitesimal surface dx a dy b . Therefore we 
can write 

f |E| = / \B\ = f |E o6 | dx a dy b = V2x Area{S) (10) 
J S J S J s 

The quantity J aa bc is the related to the angle 6 aaoc between the surface elements dx a dy b and dx a dz c . 
In fact, if we take the scalar product of the normals of these two surface elements (in the 3-space they 
span), we obtain (without writing the infinitesimals vectors) 

A ab A ac COs9aa be = 9^ {e eg hS^ ) (e fgh ,8 9 a 5^ ) = g a a9bc - 9ab9ac = Jaabc- (11) 

Finally, the 4-form 

V = ^. euKL^ 1 - 1 A E KL = ^e UKL B IJ A B KL (12) 

is easily seen to be (proportional to) the volume element dV = ^fg d 4 x. Intuitively, describing the 
geometry in terms of E rather than g ao is using as elementary variable areas and angles rather than 
length and angles. 

Using the Plebanski field, the action can be written in the BF-like form 

S[e,u>] = \\ euKL E /J [e] A F KL [u] = J B LJ [e] A F IJ [u)]. (13) 

The reason this action defines GR and not BF theory is that the independent variable to vary is the 
tetrad e, not the two-form B. While the BF field equations are obtained by varying the action (fT5|) 
under arbitrary variations of B (and uj), GR is defined by varying this action under the variations 
that respect the form ([6]) of the field B. This condition can be expressed as a constraint equation for 
E: 

S /J A jOCL = y e IJKL (w) 

Equivalently , 

*E a b ■ E cf j = 2V e a bcd- (15) 
where V = ^jVe a bcddx a A ... A dx d . This system of constraint can be decomposed in three parts: 

a) *E a6 ■ E afc = , (16) 

b) *E ah -E ac - 0, (17) 

c) *E afc -E crf = ±2V . (18) 

where the indices abed are all different, and the sign in the last equation is determined by the sign 
of their permutation. Equivalently, the B field satisfies these same equations. These three constraint 
play an important role in the following. They are called the simplicity constraints. GR can be written 
as an 50(4) BF theory whose B field satisfies the simplicity constraints (|16H17lfT8|) . 
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2.2 Selfdual structure of SO (4) 



In this section we recall some elementary facts about 50(4) and we make an observation about its 
representations that plays a role in the following. 

The group 50(4) is locally isomorphic to the product of two subgroups, each loc. isomorphic 
to 577(2): 50(4) - (577(2)+ x 577(2)_) /Z 2 . That is, we can write each U G 5*0(4) in the form 
U = (g+, g ) where g + G SU(2) + and G 5f7(2)_ and f7{7' = (g + g' + , g_g'_). This is clearly seen 
looking at its algebra 50(4), which is the linear sum of two commuting 511(2) algebras. Explicitly, let 
J IJ be the generators of 50(4). Define the selfdual and anti-selfdual generators J± := *J ± J, that 
satisfy J± = ± * J±. Then it is immediate to see that [J+, J_] = 0. The J + span a three dimensional 
subalgebra su(2) + of 50(4), and the J_ span a three dimensional subalgebra su(2)_ of 50(4), both 
isomorphic to 511(2). 

It is convenient to choose a basis in su(2) + and in su(z)_. For this, choose a unit norm vector n in 
i? 4 , and three other vectors Mj, i = 1, 2, 3 forming, together with n, an orthonormal basis, for instance 
v( = 5*, and define 

J i ± = \{*J±J)uv!n J (19) 

The 511(2) structure is then easy to see, since [J+, J+] = e l jkJ±- In particular, we can choose n = 
(0, 0, 0, 1), and v{ — Sf, and we have 

J i ± = -\e i jk J*±\j i0 . (20) 

Notice that in choosing this basis we have broken 5*0(4) invariance. In fact, the split 50(4) = 
su(2) + ©su(2)_ is canonical, but there is no canonical isomorphism between su(2) + and 511(2) _ or 
between SU(2) + and SU(2)-. One such isomorphism G n : SU(2) + — * 5t/(2)_ is picked up, for 
instance, by choosing the vector n. It sends g+ to g_, where the element (g+,g-) of 50(4) leaves 
n invariant. This isomorphism defines a notion of diagonal elements of the 50(4) algebra: the ones 
of the form aiJ l + + a^Jl. Exponentiating these, we get the diagonal elements of the 5*0(4) group, 
which we can write as U = (g,g). These diagonal elements form an S77(2) subgroup of 5*0(4), which 
is not canonical: it depends on n. It is the subgroup of 50(4) that leaves the vector n invariant. 
If we consider the 3d surface ("space") orthogonal to n, the diagonal part of 50(4) is (the double 
covering of) the 50(3) group of the ("spatial") rotations of this space; we denote it 50(3)„ C 50(4). 
Borrowing from the Lorentzian terminology, we can call "boost" a change of n. Its effect is to rotate 
the ± bases relative to one another. 

Notice that for any two-index quantity B IJ , 

-B ■ B = B\B\ + BtBt, (21) 

while 

-*B-B = B\B\-B'_B i _. (22) 

This split is independent from n, as the norms are not affected by a rotation of the basis. In particular, 
C = j J ■ J and C = j*J-J are the scalar and pseudo-scalar quadratic Casimirs of 50(4). They 
are, respectively, the sum and the difference of the quadratic Casimirs of su(2) + and su(2)_. The 
representations of the universal cover of 50(4), the group Spin{A) ~ SU{2) x SU{2) are labelled by 
two half integers The representations of 50(4) form the subset of these for which j + + j_ is 

integer. The representations satisfying j+ = j_, which clearly belong to this subset, are called simple: 
they play a major role in the BC theory as well as in the quantization below. 

The following observation plays a major role in section 14.21 Consider a simple representation 
This is also a representation of the subgroup 50(3)„ C 50(4), but a reducible one. Clearly, it 
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transforms in the representation j ® j, where j indicates the usual spin j representation of SU(2). If 
we decompose it into irreducible representations of <SO(3)„, we have 

U,j) -> j®j = 0ffilffi...ffi(2j-l)ffi2j. (23) 

The value of C4 in this representation is 2(j(j + 1)). Consider the value that the Casimir C3 of the 
subgroup SO(3) n takes on the lowest and highest-spin representations. C3 vanishes on the spin-0 
representation. On the spin-2j representation, it has the value C3 = 2j(2j + 1), which is related to 
C 4 by 

y/C 3 + 1/4 - y/2d + 1 + 1/2 = (24) 

and in the large j limit by 

C 3 = 2C 4 , (25) 

that is, if we have chosen n = (0, 0, 0, 1), 

J IJ Ju = J y J y +2J°V 0j , (26) 

which implies J 0t — 0. Therefore, the spin-zero and the spin-2j components of the simple 5*0(4) 
representation (j, j) are characterized respectively by 

spin : J i] = 0, 
spin 2j : J i0 = (27) 

in the "classical" large- j limit. 



3 Regge discretization 

We now approximate euclidean GR by means of a discrete lattice theory. A very natural way of 
doing so is Regge calculus [23]. The idea of Regge calculus is the following. The object described by 
euclidean GR is a Riemannian manifold (M, cIm), where M is a differential manifold and du its metric. 
A Riemannian manifold can be approximated by means of a piecewise flat manifold (A, g?a)j formed 
by flat (metric) simplices (triangles in 2d, tetrahedra in 3d, 4-simplices in 4d...) glued together in such 
a way that the geometry of their shared boundaries matches. Here A is the abstract triangulation and 
gJa is its metric, which is determined by the size of the individual simplices. For instance, a curved 
2d surface can be approximated by a surface obtained by gluing together flat triangles along their 
sides: curvature is then concentrated on the points where triangles meet, possibly forming "the top 
of a hill" . With a sufficient number N of simplices, we can approximate sufficiently well any given 
(compact) Riemannian manifold (M, d^), with a Regge triangulation (A,dA)0 

If we fix the abstract triangulation A and we vary d/\, namely the size of the individual n- 
simplices, then we can approximate to a certain degree a subset of GR fields. Therefore by fixing A 
we capture a subspace of the full set of all possible gravitational fields. Thus, over a fixed A we can 
define an approximation of GR, in a manner analogous to the way a given Wilson lattice defines an 
approximation to Yang-Mills field theory, or the approximation of a partial differential equation with 
finite-differences defines a discretization of the equation. Therefore the Regge theory over a fixed A 
defines a cut-off version of GR. 



1 For instance, in the sense that the two can be mapped into each other, P : M -—> A, in such a way that the 
difference between the distances between any two points, d,M(x,y) — d/\(P(x), P(y)), can be made arbitrary small with 
N sufficiently large. 
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It is important to notice, however, that the Regge cut-off is neither ultraviolet nor infrared. This is 
sharp contrast with the case of lattice QCD. In lattice QCD, the number N of elementary cells of the 
lattice defines an infrared cut-off: long wavelength degrees of freedom are recovered by increasing N. 
On the other hand, the physical size a of the individual cells enters the action of the theory, and short 
wavelength degrees of freedom are recovered in lattice QCD by decreasing a. Hence a is ultraviolet 
cut-off. In Regge GR, on the contrary, there is no fixed background size of the cells that enters the 
action. A fixed A can carry both a very large or a very small geometry. The cut-off implemented by 
A is therefore of a different nature than the one of lattice QCD. It is not difficult to see that it is 
a cut-off in the ratio between the smallest allowed wavelength and the overall size of the spacetime 
region considered. Thus, fixing A is equivalent to cutting-off the degrees of freedom of GR that have 
much smaller wavelength than the arbitrary size L of the region one considers. Since (as we shall see, 
and as implied by LQG) the quantum theory has no degrees of freedom below the Planck scale, it 
follows that a Regge approximation is good for L small, and it is a low-energy approximation for L 
largeH 

Consider a 4d triangulation. This is formed by oriented 4-simplices, tetrahedra, triangles, segments 
and points. We call v, t and / respectively the 4-simplices, the tetrahedra and the triangles of the 
triangulation decomposition. The choice of the letters v and / is dictated by the fact that in the 
dual-complex, to which we will later shift, 4-simplices are dual to vertices, and triangles are dual 
to /aceso The metric, we recall, is fiat within each 4-simplex v. All the tetrahedra, triangles and 
segments are flat (and, respectively, straight). The geometry induced on a given tetrahedron from 
the geometry of the two adjacent 4-simplices is the same. This structure will allow us to take the 
geometry of the tetrahedra as the fundamental dynamical object, and interpret the constraints implied 
by the fact that five tetrahedra fit into a single four-simplex as the expression of the dynamics. It is 
this peculiar perspective that makes the construction below possible. 

In d dimensions, a d — 2 simplex is surrounded by a cyclic sequence of d-simplices, separated by 
the d — 1 simplices that meet at the d — 2 simplex. This cyclic sequence is called the link of the d — 2 
simplex. For instance, in dimension 2, a point is surrounded by a link of triangles, separated by the 
edges that meet at the point; in dimension 3, it is an edge which is surrounded by a link of tetrahedra, 
separated by the triangles that meet at the edge; in dimension 4, which is the case that concerns us, 
a triangle / is surrounded by a link of 4-simplices Vi, v n , separated by the tetrahedra that meet at 
the triangle /. 

In Regge calculus, curvature is concentrated on the d — 2 simplices. In dimension 4, curvature is 
therefore concentrated on the triangles /. It is generated by the fact that the sum of the dihedral 
angles of the 4-simplices in the link around the triangle may be different from 2tt. We can always 
choose Cartesian coordinates covering one 4-simplex, or two adjacent 4-simplices; but in general there 
are no continuous Cartesian coordinates covering the region formed by all the 4-simplices in the link 
around a triangle. 

The variables used by Regge to describe the geometry c?a of the triangulation A are given by the 
set of the lengths of all the segments of the triangulation. Here we make a different choice of variables, 
which matches more closely what happens in the spinfoam formalism and in loop gravity. 

Consider one tetrahedron t in the triangulation. Choose a Cartesian coordinate system x a covering 



2 Since the expansion parameter A used in group field theory | 3D | is equivalent to the number of cells in Regge calculus, 
this discussion clarifies also the physical meaning of the group field theory A expansion. 

3 For coherence, we should also call e the tetrahedra, which are dual to the edges of the dual complex; this is indeed 
the common convention [3]. But in the present context this would generate confusion with the notation for the tetrad. 
We will shift to the e notation for edges only later on, in the quantum theory part. 
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Figure 1: The link of the face /, in grey, and the two group elements associated to a couple of 
tetrahedra t and t' in the link. 

the tetrahedron. Choose an orthonormal basis 

e(t) = e{{t) Vi dx a (28) 

for each such tetrahedron t. Here vi is a basis in i? 4 , chosen once and for all. This quantity of course 
transforms covariantly under a change coordinate system x a , therefore it is intrinsically defined as a 
one-form with values in R , associated to the tetrahedron. This will be our first variable, giving a 
discretized approximation of the gravitational field. 

Consider the five tetrahedra t A ,A = 1, ...5 that bound a single simplex v. The five variables 
e(tyi) are not independent, because the simplex is flat. Since the simplex is flat, we can always 
choose a common Cartesian coordinate system x a for the entire v. Let e(v) = e^(w) vi dx a denote an 
orthonormal basis describing the geometry in v. Then each e(i^) must be related to e(v) by an 50(4) 
rotation in R 4 . That is there must exist five 50(4) matrices V v t A such that 

e(v)l = (V vtA Yj e J a (t A ) (29) 

in the common coordinate patch. One can also define the transformation from one tetrahedron to the 
other: 

e(t A ) = U tAtB {v)e(t B ) (30) 
The later are of course constrained by the former: 

U tAtB (v) = V tAV V vtB (31) 

where Vt AV = V^ A . Now consider two tetrahedra t and t' sharing a face but not necessarily at the same 
four simplex and define as before the 50(4) transformation between the two Utt'- Remember that 
one can choose the same coordinate system for an open chain of four simplices linking two tetrahedra 
around the face they share. This constrains Utt' to be of the form: 

U w = V tVl ... V Vnt> (32) 

where v± ... v n is the (open) chain of 4-simplices between t and t' around their common face. Now 
one must be cautious. For two arbitrary tetrahedra sharing a face in the interior of the discretization, 
there are two such chains (see Figure [T]) To resolve this ambiguity, we make use of the orientation 
of the face /: this orientation gives a notion of clockwise and counterclockwise around the link. We 
choose the convention that Utt' denote the holonomy around the chain in the clockwise direction, 
starting at t' and ending at t. (For / on the boundary, or when one is considering a single 4-simplcx, 
of course the orientation of / need not be used to resolve this ambiguity.) 
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The arbitrariness in the choice of the orthonormal basis at each tetrahedron is reflected in the 
local 5*0(4) gauge transformations 



e(t) -» A(i) e(t), (33) 

e{v) -» A(v) e(v), (34) 

V vt -> A(v) Kt A(t)" 1 , (35) 

C/ tt , - A(t) tf w A(t')-\ (36) 

where A(t), A(u) <E 50(4) are the gauge parameters associated with the choice of basis in tetrahedra 
and vertices respectively. 

In each tetrahedron t, consider the bivector two-form 

E(i) = e(t) Ae(i). (37) 

and its dual 

B(i) - *E(t). (38) 
In components, this is B(t) = \B I J ) (t)vjvjdx a A dx b where 

B I J(t) = e IJ KL e«(t)e{;(t). (39) 

Now, consider a tetrahedron t and a triangle / in its boundary. Associate a bivector E/(i) to the 
triangle /, as follows. E/(i) is defined as the surface integral of the two- form T,(t) over the triangle 

S/(t) = ^E(t). (40) 

This is the geometrical bivector naturally associated to the (oriented) triangle, in the frame asso- 
ciated to the tetrahedron t. Its dual bivector is Bf(t) = *T,f(t). A single triangle / belongs 
to (the boundary of) the several tetrahedra t\,t2,ts,... that are around its link. The bivectors 
'Sy (ti),Y,y (t2),^j J (ts), ... are different, because they represent the triangle in the internal frames 
associated to distinct tetrahedra. But they are of course related by 

E/ J (fi) - U^JkM U tlt2 \(v 12 ) Ef L (t 2 ). (41) 

Notice that if the tetrad is Euclidean e*(f) = 5^ and the triangle is in the (1,2) plane, the only 
non-vanishing components of Y,y (t) are E^ 2 (t) = — Ey 1 (t), while the only non-vanishing components 
of Bj J (t) are B^°(t) = — £?" 3 (t). More in general, if the normal of the tetrahedron is n 1 , then for all 
the faces / of the tetrahedron t we have 

eijKL n J Bf L {t) = 0, (42) 

m T,y{t) = 0. (43) 
In particular, if we choose a gauge where n 1 = (0, 0, 0, 1), we have 

Bj?'(i) = 0, (44) 
£j?(t) = 0. (45) 

The quantities Ej J (i) and V vt can be taken as a discretization of the continuous gravitational fields 
E /J and u IJ that define GR in the Plebanski formalism. 
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3.1 Constraints on £ 



Since our aim is to promote Sj J (i) to an independent variable, let us study the constraints it satisfies. 
It is immediate to see that the four bivectors Sy^(t), £/^(i) associated to the four faces of a single 
tetrahedron satisfy the closure relation 

+ E£(t) + E£ (t) + E# (t) = 0. (46) 

Since the two- form £ /J (i) is defined in terms of a tetrad field, it satisfies the relation (|15|). or, 
equivalently, (|16ll7ll8p . Multiplying these relations by the coordinate bivectors representing triangles, 
gives the following relations. For each triangle /, we have 

* Ey(t) •£/(*) = 0. (47) 
For each couple of adjacent triangles /, /', we have 

*£/(«)•£/,(«) = (). (48) 
We call these two constraints the diagonal and off-diagonal simplicity constraints, respectively. 

The last constraint, for / and /' in the same four simplex sharing just a point, is a little subtler. 
It is 

*£/(?;)•£/,(«) = ±12V(v). (49) 

where V(y) is the four volume of the simplex and the sign depends on relative face orientation [21] . Of 
course the volume in the last formula is irrelevant. What this equation tells us is that the volume is 
the same when computed with different choices of pairs of faces in the same four simplex. Now, using 
the transformation law for the bivectors (|41| one can write it with tetrads defined on tetrahedra: 

* E/(t) • (U tt ,(v)E f ,(t')U-Hv)) = ±UV(v). (50) 

One can show that if (f3Tj) is satisfied and the first set of constraints (|46l47l48p is satisfied then (|49]) 
or, equivalently, ([50]) is satisfied automatically. 

Let's count the degrees of freedom for each four simplex, as a check. We start with the 60 degrees 
of freedom of the bivectors, 6 for each face. The constraint (|46l) imposes 24 independent equations. 
The number of independent equations of the type (|47[) are 10, and finally the constraints of the type 
(|48|) contribute 10 independent equations. The last constraint is implicitly imposed when we consider 
the tetrahedra to belong to the same four simplex. We are left with the 60 — (24 + 10 + 10) = 16 
degrees of freedom of the tetrad e(v). 

There is however some additional discrete degeneracy. The set of constraints (|46I47I48[) has two 
classes of solutions (see appendix |B|) : 

£^ 7 = 2eze£ and cyclically. (51) 

and 

E^ 7 = e IK KL e£e% and cyclically. (52) 

These two sectors of Plebanski theory are well known in the literature (see [5], [21]); they correspond 
to the fact, remarked above, that both £ and B satisfy the same equations (|46l47l48p . Because of the 
double solution ([5T]) and (fij2")) . these constraints do not determine £ uniquely. 

However, we can give the off-diagonal simplicity constraints a different and slightly stronger form, 
which fixes this degeneracy, and which will play a role below. We can replace the off-diagonal simplicity 
constraints with the following requirement: that for each tetrahedron t there exist a covariant vector 
Uj such that (|4"3"| holds for all the faces / of the tetrahedron t. It is immediate to see that this implies 
the off-diagonal simplicity constraints, and that it is implied by the physical solution (|5ip of these 
constraints. Equivalently, we require that there is a gauge in which (|45p holds. The geometry of this 
requirement is transparent: ri\ is the normal to the tetrahedron t in four dimensions, and (|43p . or 
(|45p . require that the faces of the tetrahedron are all confined to the 3d hyperplane normal to nj. 
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3.2 Relation with geometry 



How do we read out the geometry of the riemannian manifold, from the variables defined? First, it is 
easy to see that the area of a triangle / is, up to a constant factor, the norm of its associated bivector: 

y/2A f = \H f (t)\ (53) 

Notice that the l.h.s. is independent of t. This is consistent because the relation between Y,y(t) and 
T,y (t') is an 5*0(4) transformation, whence the norm is invariant, |E/(i)| = |£/(t')|. This shows that 
the definitions chosen are consistent with a characteristic requirement of a Regge triangulation: the 
boundaries of the flat 4-simplices match, and in particular the area of a triangle computed from any 
side is the same. Of course, it follows that the same is true for the volume of each tetrahedron. 

Given two triangles / and /' in the same tetrahedron t, there is a dihedral angle Off, between the 
two. This angle can be obtained from the product of the normals 

Jff ■■= AjAf, cos*,,, = E,(i) • E//(t)/2. (54) 

These are also SO (4) invariant, hence well defined independently from the tetrahedron. The two 
gauge invariant quantities Af and Jff characterize the geometry entirely. Notice that we can view 
the area as the "diagonal" part of J, ,< and write Aj = Jff. 

It is important to notice that, as mentioned in the previous section, the area is also (the square 
root of) the norm of its associated self dual, or, equivalently, antiselfdual bivector: 

A f = 2^+E)(i)+E)(i) = 2 V /-E)(i)-E)(i). (55) 
These equalities are assured by the simplicity constraints on Sj. Similarly, 

Jff, = 4+E)(t)+E),(i) = 4-S>(t)-E>,(t), (56) 

Again, these equalities follow from the simplicity constraints. 

There exist a number of relations among the quantities Jff within a single 4-simplex@ Only 10 of 
these quantities can be independent, because the geometry of a 4-simplex is determined by 10 numbers. 
In particular, all angles must be given functions Jff(Af) of the 10 areas (up to degeneracies). Explicit 
knowledge of the form of these functions would be quite useful in quantum gravity. 

Finally, consider now a triangle /. Let t\,t2, ■■■,t n be the set of the tetrahedra in the link around 
/ and «i2, «23, •••jKni be the corresponding set of simplices in this link, where t% bounds vn and V23 
and so on cyclically. In general, if we parallel transport e(t) across simplices around a triangle, using 
Utt' , we come back rotated, because of the curvature at the triangle (the analog of a parallel transport 
around the tip of a pyramid in 2d.) In other words, we can always gauge transform e^(t) to 5^ within 
a single cartesian coordinate patch, but in general there is no cartesian coordinate patch around a 
face. We define 

C7>(ti) = V tlVl2 V Vl2t2 ... V tnVnl V Vnltl (57) 

or, equivalently 

U f (h) ee U tlt2 (v 12 ) ... U tntl (v nl ), (58) 

the product of the rotation matrices obtained turning along the link of the triangle /, beginning with 
the tetrahedron t\. The rotation matrix Uf(t\) then represents the SO(4) curvature associated to the 
triangle /, written in the ii-frame. 



4 As a discussion of these relations does not seem to appear elsewhere in the literature, we have here included an 
appendix on them JaJ. 
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3.3 Dynamics 



The last step before attacking the quantization of the model is to write the discretized action. Take 
the e(t) as independent variables so that £/(i) = e(t) Ae(i). In analogy with Regge calculus we define 
the action to be 

S bu i k [e(t),U,V] =\Y. Tr \ B f( t ) u f( t )]+Y.Y. Tr ^f Utt '^ Vt 'v ( 59 ) 

/ v /Co 

where we recall that B = *S. The first sum is over all faces in the interior of the discretization and 
Uf(t) is defined as in (|58p, where t is any one of the tetrahedra in the link of /. We impose a priori 
equation (|4Tj) , which implies @ 

Xf(t)U tt > =U a >X f (t'). (61) 

It follows that the first term in the action |59|) is independent of the choice of tetrahedron t in the 
link of each /. 

The second term in (|59[) is a sum over all four simplices and in each four simplex over all faces 
belonging to it@. X v f is a lagrange multiplier living in the algebra of SO (4) and varying with respect 
to it gives the constraint ([3"Tj) on the group variables. 

This action is invariant under the gauge transformations (|36p and reduces to the action of GR 
in the limit in which the triangulation is fine. This can be seen as as follows. In the limit in which 
curvatures are small, Uf(t) = 1 + ^F a t>dx a A dx b , where dx a A dx h is the plane normal to the triangle 
/. Hence the trace gives 

l -Tr[B f {t)U f {t)] ~ \Bl J b F c K d L e abcd = \e IJ KL e* e L b F* L e abcd = ee&Ftf = JgR, (62) 

which is the GR lagrangian density. 

There is also a close relation of ([59]) to the Regge action. To see this, let us first see how to extract 
the deficit angle around each face in our framework. Let v±, v 1 ^ denote two of the edges of the triangle 
/. Now, as the triangle / is the axis of the parallel transport around /, the tangent space parallel 
transport map Uffy^v preserves each of the vectors v^v^: 

£Mf)^< = <^e{l,2}. (63) 

Contracting both sides with e*, and inserting the resolution of unity e u j{t)e p (t) = 8" we get 

U f (tYjvi =v I i , te {1,2} (64) 

where v\ :— e^v^ . By linearity, Uf(t) I j therefore acts as identity on the subspace V := span{i>(, 
forcing C/y(t) 7 j to belong to the SO(2) subgroup fixing V. Therefore Uf(tY j is described by a single 
angle: this is the deficit angle around /. To make this explicit, as well as to cast (|59"|) directly into 



5 It follows that, for each /, there is only one independent Bf(t) in the link. This is consistent with the action l|59|l . 
where only one Sj-(t) in each link appears. Note it follows that we are also imposing a priori 

Z f (t)U f (t) = U f m f (t), (60) 

so that not even this one independent in the link is a priori independent of the connection degrees of freedom. 

6 This is equivalent to summing over the wedges defined by Reisenberger in |16| . Each wedge is identified by a pair 
(*./)■ 
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Regge-like form, let us introduce an orthonormal basis £| of the orthogonal complement of V. The 
rotation Ufit) 1 j is a rotation in the £i-£2 plane. Therefore 

U f (t) 21 := (6), UfdYjil = S mO f (65) 

where 9f is the angle by which £1 is rotated by Uf(t) — the deficit angle. We have also 

U f (t) 12 := {Ci)iU,(tYjCi =-sin9 f . (66) 

Furthermore, 

Bf(t) 1J ocgtf. (67) 
The norm of Bf(t) IJ is just the area (|53|) . which fixes the constant of proportionality here. We have 

B f (t) IJ =A f gg. (68) 

Thus 

B f {t) IJ U f {t)u = A f g$Uf(t)ij = 2A f sme f . (69) 
The bulk term in (|59|) can therefore be written 

^2 Af sin Of . (70) 
/ 



To lowest order in 6t, this is precisely the bulk Regge action. 



3.4 Classical equations of motion 



From the form of the action ([62]) . one can see that variation w.r.t. the tetrads will give the discrete 
analog of the Einstein equations. Variation with respect to the connection is more subtle. In order 
to proceed, let us write the action on shell w.r.t. the constraints on the group variables so that it 
reads S[e, V] = ^J2f Tr [ B f( t ) u f( t )} where U f (t) is written as in (57|). The action, as we saw, is 
independent of the choice of the base tetrahedron in each face. Let us consider the variation with 
respect to Vt v ■ This variable appears in four terms in the sum, corresponding to the four faces of the 
tetrahedron t. Let us choose in addition, for the sake of simplicity, the tetrahedron t to be the base 
of these four faces. The variation (SVtv = £Vt„)0 gives 

6S= l E Tr[B fi (t)ZU fi (t)} (71) 

/i,i=1..4 

where £ is an arbitrary infinitesimal clement in the algebra of 5*0(4). Stationarity w.r.t. these 
variations implies that the antisymmetric part of Ylf- ^ h^h > seen as a f° ur dimensional matrix, is 
zero. Explicitly, 

Y,U U B U +B h Uj l =Q. (72) 

ft 



7 The variation so defined is that determined by the right invariant vector field associated to £. 
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Now, U f (t) can be expanded as Uf ~ 1 + Uf where Uf is of second order in the lattice spacing 0, 
so that to first order we get simply, 

which is the closure constraint. The analogous continuum equation is the Gauss constraint, given by: 



DB = dB+ [A, B] = (74) 

where A is the spin connection. Remembering that the tetrahedron is flat, one can choose the con- 
nection to be identically zero. Integration over the region defined by t gives the closure. In order to 
make this relation more clear, let us consider now the case where t is not the base tetrahedron for all 
faces. Consider then, for the face /i the base tetrahedron to be t\, where t,t\, t n is the link around 
/i. Variation w.r.t. V tv gives: 

[U m B h (h)U^U fl (t) + ... + B h (t)U u ] ant . = (75) 
where Utt x — Vt Vl V Vl t 1 - To first order it reads: 



B fl (h) + [uu^Bf.ih)} + ... + B h (t) = (76) 

which can be seen directly as the integration of the equation (1741) over the region defined by t, where 
the spin connection is a distribution concentrated at the face f\ in the interior of t. 



3.5 Topological term 

Recall that in the LQG approach the action that is quantized is the Hoist action [33] , obtained adding 
to the action a topological term that doesn't change the equations of motion 

S Ho ut= J *(eAe)AF + - J (eAe)AF. (77) 

where 7 is called the Immirzi parameter. The introduction of a topological term is required in order 
to have a theory of connections on the boundary: without it, as shown by Ashtekar, the connection 
variable does not survive the Legendre transform [34] . It is easy to add a similar term in the discrete 
theory, having no effect on the equations of motion in the continuous limit. This is 

S top [e(t), U, V] = i E Tr[V f {t)Uf{t)]. (78) 
7 / 

The reason this has no effects on the equations of motion is interesting. In a Regge geometry, the 
curvature associated to the face / is given by the rotation Uf(t). This rotation has the property of 
leaving the face / itself invariant. Hence, it is a rotation that is generated by the dual of the face 
bivector £/(£) itself. That is, it has the form Uf(t) ~ exp{9Bf(t)}. In the weak field limit, this gives 
U f (t) - 1 + 6B f (t), and therefore Tr[E f (t)U f (t)] - 6Tr[T,f(t)B f (t) = 6Ef(t)*Ef(t)] which vanishes 
because of the simplicity constraint (because £ is simple). 

This term will play a role in the quantization. 



8 Latticc spacing is defined with respect to a particular continuum limit. To be specific, suppose one has a metric g 
which one wishes to approximate. Then, for each real number A, one can construct a corresponding Regge geometry 
approximating g, such that the typical lattice spacing is A, as measured w.r.t. g. As A — » 0, approaches g, and 
one can consider expansions in A. 
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Figure 2: The face / and the two tetrahedra t and t 1 are on the boundary of the triangulation, 
represented by the heavy dashed line. The link of the face is divided by the boundary into two parts. 

3.6 Boundary terms 

In ordinary Regge calculus boundary terms must be added to the action so that the equations of motion 
are the same when we vary the action with some boundary variables fixed. The other condition for 
boundary terms is that they add correctly. Consider a region of spacetime which is the union of 
two disjoint regions separated by a boundary, S = S' U S" . The additivity condition is then that 
S[S] = S[S'] + S[S"] ps]. 

Consider for simplicity the part of the action referring to a single face Sf — Tr[B f(t)U /(£)]■ 
Furthermore, choose two tetrahedra t and t' appearing in the link around this face and break the link 
in two so that the two tetrahedra belong now to the boundary (cf. Figure [5]). The action can be split 
to first order in the algebra as 

Sf = \Tr[B f {t)Ul t ,R f ] + ^TrlBfit'^Rj 1 } (79) 

where in the second term of the sum we have replaced Bf(t) by Bf(t') which, because of (|4lj). 

(12). 

introduces terms of the second order in the lattice spacing. U^/ are just the boundary connection 
variables. Here Rf is a fixed SO(4) element inserted to make the terms SO(4) covariant. Thus, the 
boundary action can be written in general as 

S B a = IY, Tr[B f {t)U w R f ] (80) 

/C9A 

with U U ' given by 

U tt > =U tt M) - U tnt '(v n ). (81) 
This gives additivity to first order in the algebra. 

Let us remark on the specific case when the -B's are fixed on the boundary. In this case, in the 
continuum Plebanski theory, the boundary term is just J B A A where A is the connection 36J. If we 
choose Rf = 1, the above discrete boundary action becomes 

S ^ = \ E Tr[B f {t)U w ] (82) 

fCdA 

which reduces precisely to the continuum boundary action. Thus we choose Rf = 1 when fixing 
the B's on the boundary. The boundary action is not gauge- invariant, as it is not in the continuum 
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theory^ 

Let us now write the action for a single four simplex as it will be useful to us in the next section. 
All the faces are on boundary and so the action is just a sum of boundary terms. Explicitly: 

S " = \ E Tr i B f(tA) U tAtB (v)} + TV [A/ U tAtB (v) V tBV V vtA ] (83) 

fCv 

3.7 Boundary variables 

Suppose the triangulation A has a boundary dA. This boundary is a 3d manifold, triangulated by 
tetrahedra separated by triangles. Notice that, unlike what happens in the bulk, each boundary 
tetrahedron bounds just a single simplex of the triangulation; and each boundary triangle bounds just 
two boundary tetrahedra. 

Let us identify the boundary variables, which reduce to the boundary gravitational field in the 
continuum limit. One boundary variable is simply £/(t), where / is a boundary triangle and t is 
a boundary tetrahedron. There are only two boundary tetrahedra around /, one, t, at the start of 
the link as determined by the orientation of /, and the other, i', at the end. It will turn out to be 
convenient to use the notation and E^ for E/(i) and E/(i'). If desired, one can also associate 
variables with the reverse orientation of /: ^f-i '■= Ey , ^j-i '■= Ey . 

The other boundary variable is the group element Uf = U t t> giving the parallel transport across 
each triangle / bounding t and t 1 (not to be confused with the holonomy ([57)) defined above) . Notice 
that (HJ) implies 

U f Ef Uj 1 = E£. (84) 

Finally, of the constraints flU E3 03 S3 ED) , only ([H E3 g3 gSJ act separately at each tetrahe- 
dron, and thus impose direct restrictions on boundary data. We call ([41)1 14"3 14"51 14"3"| the kinematical 
constraints. The constraint ([3Tj) on the other hand necessarily involves bulk variables; we thus call it 
the dynamical constraint. 

The complex dual to the boundary triangulation defines a 4-valent graph T with nodes t (dual to 
the boundary tetrahedron t) and links / (dual to the boundary triangle /). It is convenient to view 
the boundary variables as associated to the four valent graph T dual to the boundary triangulation: 
we have SO(A) group elements associated to the links of the graph and two E variables associated 
to the two orientations of each such link. Accordingly, we change notation, and call I (for links) the 
oriented boundary triangles and n (for nodes) the boundary tetrahedra. 

Since we are in a first order formalism, the space of these variables (T,f,T,f,Ui) code the phase 
space of discretized GR. In fact, this space is precisely the same as the phase space of a Yang-Mills 
SO (A) lattice theory, and it can be identified as the cotangent bundle of the configuration space 
C = SO(A) L , where L is the number of links on T. 

This cotangent bundle has a natural symplectic structure, which defines Poisson brackets. 

{U U U V } = 

{(Ef) /J ,£M = fo, tfjr" 

{(Ef) /J ,£M = Sw^Ui 

{(Ef)",(Ef)^} = 8 w \ IJKL m N {^) mn (85) 



9 One can also fix the U's on the boundary. In this case, there is no boundary term in the continuum theory. This 
can effectively be achieved in the discrete theory by setting Rf = U~ t } after varying the action to obtain the equations 
of motion (see appendix [Cj. 
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where t ij and X\j^ L are, respectively, the generators and the structure constants of 5*0(4). Observe 
that the two quantities E^ and E Z L act very nicely as the right, and respectively, left invariant vector 
fields on the group. Equation ([84")) gives the correct transformation law between the two. 

However, it is important to notice a crucial detail at this point. As pointed out in [3§] and in [25] , 
because of its peculiar SU(2) x SU{2) structure, the cotangent bundle over 5*0(4) carries indeed two 
different natural symplectic structures, related to one another by what Baez and Barrett call a flip: 
one is obtained from the other by flipping the sign of the antiselfdual part. That is, replacing E by B 

{U h U v } = 

{{Bl) IJ M = S w U lT IJ 

{{B^) IJ ,U V } = S W T IJ Ut 

{(B*) IJ ,{B*) KL ] = 8 W \ IJKL MN {B*) MN (86) 

Both structures are equivalent to the lattice Yang Mills theory Poisson brackets [32] , the difference is 
only whether we identify the electric field with B or with E. We call (|85|) the "flipped" Poisson struc- 
ture (at the risk of some confusion, since Baez and Barrett call (|8"6")l the "flipped" Poisson structure. 
Being flipped as everything else, is a relative notion...) 

Which one is the correct Poisson structure to utilize? The classical equations of motion, which 
is the part of the theory that is empirically supported, do not determine the symplectic structure 
uniquely. In the Appendix [C] we study the direct construction of the symplectic structure from the 
action. If we take the action ([55]) without the topological term ([75]) . then we have the unflippcd 
Poisson structure ([86]) . (This is why we call it "unflipped" .) But in order to arrive at the Ashtckar 
formalism and LQG, we know that the topological term is needed. As shown in the appendix, the 
two symplectic structures written above are recovered by taking 7> 1 and 7 -C 1, respectively. We 
discard the first choice that gives macroscopic discrete area eigenvalue and we choose the secondf^l 

Thus, we choose the flipped symplectic structure (|85|) as a basis for the quantization of theory. This 
choice leads to a nontrivial intertwiner space, while the opposite choice leads to the Barrett-Crane 
trivial intertwiner space, as we shall see below. 

3.8 Summary of the classical theory 

Summarizing, discretized GR can be defined by the action (|59| with the appropriate boundary terms 
(|80p now defined as a function of the variables [Sf(t), V v t, Utt>(v), ^ v f})- That is 

5 6uifc [E, [/,!/, A] = \J2 Tr l B f( t ) u f( t )]+J2J2 Tr ^f Utt '^ Vt 'v V «A- ( 8? ) 
/ « fCv 

Plus the diagonal and off-diagonal simplicity constraints 

C ff := \*X f (t) ■ E/(t) = 0, (88) 
C//' :=!*£/(*)•£/'(*)= 0.. (89) 

The second can be replaced by the (stronger) condition that there is an nj for each t such that 

m Zf J (t) = 0, (90) 



10 Intermediate choices may also be of interest, but they lead to a more complicated quantum theory, involving also 
non-simple representations. See Appendix Icl 
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equivalently, there is a gauge in which 

Tfp(t) = 0. (91) 
The following two other constraints follow from the equations of motion 

= 0, (92) 

and 

U tAtB (v) = V tAV V vtB . (93) 

4 Quantization 

The quantization of the theory will proceed in three steps. First, we write the Hilbert space and the 
operators that quantize the boundary variables and their Poisson algebra. Second, we impose the 
kinematical constraints (|92l88l89l9ip . The first two of these pose no problem. The third one, namely 
the nondiagonal simplicity constraints ([89]) require a careful discussion and some technical steps. As 
we shall see, we cannot impose these constraints strongly; we impose them weakly, in a sense defined 
below. Finally, the dynamics is specified by computing an amplitude for a state in the boundary 
Hilbert space [3TI [35]. This amplitude is constructed by building blocks, the elementary building 
block being a single four-simplex. 

Remarkably, we will find that the physical Hilbert space that solves the constraints is naturally 
isomorphic to the kinematical Hilbert space of SO(3) LQG. More precisely, it is isomorphic to the set 
of states of SO(3) LQG defined on the graph formed by the dual of the boundary triangulation. Our 
key result will then be a transition amplitude for a state in the Hilbert space of 50(3) spin networks 
defined on the graph corresponding to the (dual of the) boundary of a four-simplex. 



4.1 Kinematical Hilbert space and operators 

The boundary phase space is the same as the one of an 50(4) lattice Yang-Mills theory. We quantize 
it in the same manner as SO(A) lattice Yang-Mills theory [35]. That is, the natural quantization of the 
symplectic structure of the boundary variables is defined on the Hilbert space 7iso(4) := ^2[SO(4) L ], 
formed by wave functions ip(Ui) = ip{g^~ 7 S; - )- The XJ\ variables are quantized as diagonal operators. 
The variables Yjf and are quantized as the left invariant — respectively right invariant — vector 
fields on the group. These define a representation of the flipped classical Poisson algebra ([55)1 . 

A basis in this space is given by the states 

*itiTr*iM>9T) = lstaT\3i>3T>%,K) = (®DW(gt) ■ ® J+j fester) ■ (8)i» J ■ 

(94) 

Here D^\g) the are matrices of the SU(2) representation j: two of these are associated to each link; 
together, they form the SO(A) representation matrix (j + , defined on the Hilbert space T~L(j+,j-y 
At each node, Jj" € "H-+ ® • • • ® H.-+ and ir £ H-- <8> • • • ® TL-- , so that 1+ ® I~ defines an element 
of the space 

n Ut,37)-Ut,I) : ( H Ut,H) ® - ® H -u ,, ( 95 ) 

associated to each node n with fixed adjacent representations (jf , j\)---{jt , j^)- The contraction 
pattern of the indices between the representation matrices and the tensors is defined by T. 

Next, we consider the constraints. 
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Strictly speaking, the closure constraint (|9"2")1 does not need to be imposed at this stage, since it is 
not an independent constraint, but it is implemented by the dynamics, as shown in the last section. 
But, anticipating, let us see what is its effect. It reduces L2[SO(A) L ] to the space L2[SO(4) L / SO(A) N ] 
formed by the states invariant under ip{Ui) = i>{V s {t)UiV^) . That is, it constrains the tensors (/ + , I~) 

to be SO(A) intertwiners That is, it constrains I + <g> I~ to belong to the subspace /cf°^ := 

. + := Inv ( H (i + — ■) <X> ... ® H, + ,—•)), which is the space of the SO(A) intertwiners. The 



state space obtained by imposing the closure constraint is then precisely the Hubert space of an SO (4) 
lattice Yang-Mills theory. 

The diagonal simplicity constraint (|88p restricts to the spin network states where the representation 
associated to the links is simple. That is, it imposes = j[~ = jt . 

Let us now come to the off-diagonal simplicity constraints (|89|) . which are of central interest to 
us. After (|92|) and (|88|) are satisfied, only two of the three off-diagonal simplicity constraints acting 
on each tetrahedron are independent. These constraints form a second class system. Imposing them 
strongly restricts the space of intertwiners to one unique solution given by the Barrett-Crane vertexP^I 

In order to illustrate the problems that follow from imposing second class constraints strongly, 
and a possible solutions to this problem, consider a simple system that describes a single particle, but 
using twice as many variables as needed. The phase space is the doubled phase space for one particle, 
i.e., T*M. x T*R 9 Pi), (92,^2)), and the symplectic structure is the one given by the commutator 
{la,Pb\ — S a b- We set the constraints to be 

q 1 -q 2 = 0, 

Pi-J»2 = 0. (96) 

By defining the variables q± — (qi ± <?2)/2 and p± = (pi ±p2)/2, the constraints read: q_ = p_ = 0. 
They are clearly second class since {q~,p-} = 1. Suppose we quantize this system on the Schrodinger 
Hubert space L,2[R 2 ] formed by wave functions of the form ip(q+, <?-). If we impose the two constraints 
strongly we obtain the set of two equations 

Q- <?-) = 0, 
ih-^-^(q + ,q_)=0. (97) 

which has no solutions. We have lost entirely the system. 

There are several ways of dealing with second class systems. One possibility, which is employed for 
instance in the Gupta-Bleuler formalism for electromagnetism and in string theory, can be illustrated 
as follows in the context of the simple model above (see for instance the appendix of [Mj)- Define the 
creation and annihilation operators a1_ = (p_ + iq-)/y/2 and a_ = (p_ — iq-)j\[2. The constraints 
now read a_ = aL = 0. Impose only one of these strongly: a_|-0) = and call the space of states 
solving this Ti p h- Notice that the other one holds weakly, in the sense that 



■|ot|^) = V^ieV (98) 

That is, a^L maps the physical Hilbert space Tt p h into a subspace orthogonal to Tt p h- Similarly, in the 
Gupta-Bleuler formalism the Lorentz condition (which forms a second class system with the Gauss 



11 The commutator of two of these constraints is called chirality in [7]. In [7], the chirality constraint is imposed 
strongly on the states as well, with the result of selecting the non-degenerate geometries corresponding to the sector 
151> . Notice, however, that the system formed by the the chirality and the simplicity constraints is not first class either, 
as the chirality, in turn, does not commute with the simplicity constraints. 
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constraint) holds in the form 

{4>W l A^)=0 V<j>,^£H ph . (99) 

A general strategy to deal with second class constraints is therefore to search for a decomposition 
of the Hilbert space of the theory Hkin — Hphys © 7~tsp (sp. for spurious) such that the constraints 
map "Hphys ~~ > H S p. We then say that the constraints vanish weakly on TL p h ys - This is the strategy 
we employ below for the off-diagonal simplicity constraints Ciy . Since the decomposition may not be 
unique, we will have to select the one which is best physically motivated. We now define this space. 



4.2 The physical intertwiner state space JC p h 

Consider the state space obtained by imposing the diagonal simplicity constraint, namely by taking 
jj + = jT~ := ji, but not the closure constraint yet. Let us restrict attention to a single tetrahedron 
t. The constraints (|59")l at t act on the space associated to i, which is 'H(j 1 ,j 1 )...(j i j 4 ,)- The closure 
constraint will then restrict this space to the SO (4) intertwiner space /Cf . We search for a subspace 
fCph of /cf where the nondiagonal constraints vanish weakly. 

First, note K.f°^ = liw(H^j 1 j 1 ) <g> ■■■'H(j i ,j 4 )) is a subspace of the larger space 

® -^(J4,i4) = "Hh®h ® -Wj 4 ®i4 ( 10 °) 

which can be thought of as a tensor product of carrying spaces of 50(4) representations or SU{2) 
representations, as desired. The Clebsch-Gordan decomposition for the first factor on the right-hand 
side above gives 

^h9h = n h ® H h = Ko ® Hi © • • • © H 2jl (101) 

and similarly for the other factors. By selecting the highest spin term in each factor, we obtain a 
subspace 

H 2jl ®---®H2 U - (102) 

Orthogonal projection of this subspace into K,^° then gives us the desired K p h C K-f°^ . This IC p h 
is the intertwiner space that we want to consider as a solution of the constraints. The total physical 
boundary space TC p h of the theory is then obtained as the span of spin-networks in L 2 [SO(4) L / SO(A) N ] 
with simple representations on edges and with intertwiners in the spaces K, p h at each node. Notice 
that the elements in IC p h are not necessarily simple in their internal representation, in any basis. Let 
us study the properties of the space IC p h, and the reasons of its interest. 



(i) First, it is easy to see that the off-diagonal simplicity constraints (f59"| vanish weakly K p h, in the 
sense stated above. ' 
can be expanded as 



sense stated above. This follows from the following consideration. A generic element of /cf ° 



^ c i+i -\i+,r) (103) 



where defines a basis in the SU (2) intertwiner space. The off-diagonal simplicity constraints 
are odd under the exchange of i + and i~ , namely under exchange of self dual and antiselfdual 
sectors. But the states in K, p h are symmetric in i + and i~ . Hence {4>\Cw \ip) = 0, V0, if) € IC p h, 
that is, IC p h can be considered as one possible solution of the constraint equations. 

(ii) Second, let us motivate the choice of this solution. Recall that the off-diagonal simplicity con- 
straints can be expressed as the requirement that there is a direction nj such that (f9~Tj) holds. 
But this is precisely the "classical" limit of the condition satisfied by the spin-2j representation, 
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as observed at the end of Section |2~21 That is, promoting (p?Tj) to the quantum theory gives the 
requirement that there there is a gauge in which 

J 0i = 0. (104) 

which is (|27[) . namely the condition satisfied by the spin-2j component of the representation 
Equivalently, (f9l"j) implies 

2 ^4 = \j{ J Ju.j = \jfJia = C 3 (105) 

where C4 is the 50(4) scalar Casimir in the representation Ti.j l j t and O3 is the Casimir of the 
50(3) subgroup that leaves n 1 invariant, in the same representation. As it is, this relation 
has in general no solution in tLj l j l , but if we order it in slightly different manner, as in (|24|) . 
(reinstating h 7^ for clarity) 

C = VCs + h 2 /4 - \/2C 4 + K 2 + h/2 = Q (106) 

then there is always a solution, which is given by the H^j subspace of Hj.j- Therefore the off- 
diagonal simplicity constraints pick up precisely the space we have defined. In other words, this 
space satisfies a quantum constraint, which in the classical limits becomes the classical constraint 
GUI). 

Notice that if we hadn't chosen the flipped symplectic structure, then we would have obtained 
J % i = 0, instead of J° l = 0, that is, the vanishing of the 50(3) Casimir. Following the same 
procedure as above, we would have defined the space 

4° } = Inv (Ho <8 ... ® n ) C H n ... u (107) 

and this is precisely the one dimensional Barrett-Crane intertwiner space. That it, it is the 
choice of the flipped symplectic structure that allows the selection of a nontrivial intertwiner 
space. 

(iii) Third, we have the remarkable result that IC p h is isomorphic to the £0(3) intertwiner space, 
and therefore the constrained boundary space H p h is precisely the 50(3) LQG state space 
H.so(3) associated to the graph which is dual to the boundary of the triangulation, namely 
the space of the 5*0(3) spin networks on this graph. We exhibit this isomorphism in a way 
that simultaneously shows a new way of viewing TL p h- We will first construct a projection 
T : ^so(4) * ^so(3)- The hermitian conjugate map will then be an embedding / : Tiso(S) ~ * 
H.so{A)- The image of this embedding will be none other than Tl p h- 

The map n : Hso(4) — * ^so(3) is simple to construct. Choose an S*0(3) subgroup H of 50(4). 
As explained in Section 12. 2[ this choice is equivalent to the choice of an isomorphism between 
the two SU(2) factors of 5*0(4). A (simple) 50(4) representation TCtjj) is in general a reducible 
representation of H: since H acts on each duality component independently, it transforms as 
Hrj j) = Hj (8> Hj, where Hj is the standard spin-j representation of SU(2). This decomposes 
into irreducibles as 

h U:J) = n e Hi e ... e n 2j , (ios) 

Thus, begin with an 50(4) spin-network state Vf, and consider the restriction W| ; to the 

xfl 

subgroup H on each edge. As ^ is 50(4) invariant, it is in particular -ff-invariant, and thus the 
restriction \£| , is a sum of 50(3) spin- networks. Furthermore, because 4" is 50(4) invariant, 

and because all possible 50(3) subgroups H are related to each other by conjugation, this sum 
of spin-networks is independent of the choice of 50(3) subgroup H . The above considerations 
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tell us that this sum is just that given by the above Clebsch-Gordan decomposition for each edge. 
The result of the projection ir^S is then defined to be the spin-network in this sum corresponding 
to the highest weight spin on each edge. Thus tt so-defined maps simple SO (A) spin-network 
states to 50(3) spin-network states. The spin on each edge of the 50(4) spin-network is 
mapped to spin 2j in the SO (3) spin-network. Each 5*0(4) intertwiner I v e lny SO u)('Hrj lt j 1 -\ ® 
... ® ^(j 4 ,j 4 )) is mapped to an SO(3) intertwiner simply via orthogonal projection onto the 
subspace Inv sc/(2 )(W( 2 j 1 ) ® •■• ® W(2j 4 )) C (H {2jl ) ® ... ® H (2j4) ). 

Let us now describe the conjugate embedding / : Ttso(3) ~^ 7~tph- This is defined as the hermitian 
conjugate of it, using the fact that tt is a linear map between Hilbert spaces. Let us describe it 
in detail. 

Let us first describe the embedding / restricted to a single intertwiner space, namely the map 
(that we also call) / from the space IC SO (3) of the 5*0(3) intertwiners to JC p h- Consider an 
intertwiner i G A^so(3): between the four representations {2j\...2j±). Contract it with four 
trivalent intertwiners between the representations (2j a ,j a ,ja), the edge with spin 2j\ being 
contracted with the 2ji edge of the corresponding trivalent intertwiner, etc. This gives us a 
tensor e(i) in ® •■• ® ^(j 4 .j 4 ))- e (*) i s n °t an 50(4) intertwiner, because it is not 

50(4) invariant, but we obtain an 50(4) intertwiner by projecting it in the invariant part of 
{'H(j 1 j 1 ) ® ... ® 7~L(j 4 ,j 4 ))- Since 50(4) is compact, this projection can be implemented by acting 
with a group element U in each representation, and integrating over 50(4). 




(109) 



The 50(4) action can be factorized into two SU (2) group elements, one acting on the self dual, 
and other on the anti-selfdual representations. One of the two factors can be eliminated by 
virtue of the SU{2) invariance of the trivalent intertwiners and i. What remains is an SU{2) 
integration over just one of the representations. Using the well known relation 

/ dgD aibl (g)...D aibi (g) =S2i ai - a H bl - M (110) 

JSU(2) i 

it is easy to see that we have 

m^Ystt^-v+n an) 

where the coefficients are given by the evaluation of the spin network 




If we piece these maps at each node, we obtain the map / : Hso(3) ~^ 7~Lph of the entire LQG 
space into the state space of the present theory. In the spin network basis we obtain 

/ : \j u i n ) ^ J2 /^Vl^2,i//2,*+,0- (H2) 

it, in 
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Equivalently, writing explicitly the states as functions on the groups, 




l[dV n ® I**'"*) (V aW {gt>9r)V t (i) ■ ® (113) 

'SO(4) N n V i / V n J 

where indices have been omitted and s(l),t(l) stand resp. for source and target of the link I. 
This completes the definition of the projection and the corresponding embedding of the Hilbcrt 
space of LQG into the boundary Hilbert space of the model. 

) Let us illustrate more in detail the construction in (iii) using the standard spinor notation 
|40j . The vectors in Tij can be represented as totally symmetric tensors with 2j spinor indices 
(Ai...A-2j) = A, where each index A — 0, 1 is in the fundamental representation of SU(2). An 
element in 7~t(j 1 .j 1 )...(j 4 ,j 4 ) has therefore the form 

I (A 1 ...A 2jl ){A' 1 ...A' 2jl )...(D 1 ...D 2U )(D' 1 ...D' 2ji ) _. jAA'.-.VV' q 14 n 

here primed and unprimed indices are symmetrized independently; they live in the self dual and 
antiself dual components of the representation. Round brackets stand for symmetrization. By 
choosing to no longer distinguish between primed and unprimed indices, this 50(4) intertwiner 
I between the simple representations (ji, Ji)...(j4, J4) becomes a tensor among the SU(2) rep- 
resentations ji <g> ji,]2 ® J2, J3 ® hiji ® J4- Because of the SO(4)-invariance of /, the resulting 
SU(2) tensor does not depend on the way primed and unprimed indices arc identified. 

Let us first construct the projection 7r, and corresponding embedding /, for individual intertwiner 
spaces. The projection 7r : ICsoU) ~ * ^so(3) i s simply given by symmetrizing over the spinor 
indices associated with each link 

n . jAA'...VV' , y j(AA')...(DV') _. ^...d 

where the index a is short for (A/4/) := A\ . . . A2j 1 A[ . . . A' 2 j 1 , and similarly for b,c, and d. This 
symmetrization projects / to an SU(2) intertwiner among the representations 2ji, 2j2, 2j3, 2j4, 
thereby selecting the highest weight irreducible representation in the decomposition j ® j = 
® • • • 2j on each edge. Notice that 2j e Z so that the projected intertwiner transforms under 
SO(3) transformations. 

Next, we write, for ji,...j4 € Z, the corresponding embedding / : ^50(3) — > TCgoU)- First, 
the embedding e is trivially obtained by reading the indices AA' as living in self-dual and 
anti-self-dual representations, respectively: 

e : i a ~ d e(i) AA '-' DV ' := i^'-™' (116) 

This is not yet a SO(4) intertwiner and in order to recover an element of the invariant subspace 
we have to "group average" as in equation (|109jl . 

Let us now consider the projection, and corresponding embedding, for individual spins. Consider 
the state space for a single link - that is, the space of square integrable functions over one single 
copy of 5*0(4). This space is spanned by the representation matrices D^> J •* (Ui) AA bb' where 
Ui = (g?,gf) € SO (A) and gf 1 € SU(2). Here we restrict ourselves to simple representations 
For such representation matrices, the projection is defined by 

tt : DM{g+,gr) AA 'BB> — > £ (iJ) (ffi, <?0 ( ^'W) = D™{gi)%. (117) 
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That is, we first restrict to a diagonal subgroup {(g,g)} C SO(4); such a subgroup is isomorphic 
to SO(3). The SO(4)-representation matrix then becomes a j ® j SU(2)-representation 

matrix; from there we project to the highest SU(2) irreducible, namely the 2j representation, 
as before. The corresponding embedding of SO(3) spins into SO(4) spins is then obviously 

Putting together the embeddings for intertwiners and spins, we obtain the embedding (|113[) of 
5*0(3) LQG states into the 50(4) boundary state space of the model. 

(v) Lastly, with the embedding proposed in (iii,iv) above, and in light of (|53"|) . the constraints (| 1 05[) 
being used to solve the off-diagonal simplicity constraints simply express the condition that 
three-dimensional areas as determined by the SO (4) theory match three-dimensional areas as 
determined by LQG. 

This concludes the discussion on the implementation of the constraints. The last point to discuss is 
the dynamics. 

4.3 Dynamics 

Following the strategy stated at the start of this section, consider a triangulation formed by a single 
4-simplex v. Denote the boundary graph by r 5 . A generic boundary state (satisfying all kincmatical 
constraints) is a function ^>(U a b), where a, b = 1, 5, in the image /[Hlqg]- We begin by writing the 
transition amplitude between sharp values of the B's. From this we deduce the amplitude in terms of 
the U's, and then compute the amplitude for the quantum state ^(U a b). Beginning this procedure, 



A(B ab ) = 




(118) 



so that the amplitude in the connection representation reads: 

MUab) = I Y[dV a SiVaUabV^ 1 ). (119) 
J a 

The integral is over a choice of S0(4) element V at each node. Let us define the bra 




(120) 



The quantum amplitude for a given state \ip) in the boundary Hilbert space is: 

:= (W\V) (121) 

As *f> satisfies all the kinematical constraints, it is of the form $ = f(ip) for some ip £ H-lqg- Let us 
consider the specific case when ip is a spin-network state ip = ^{j a b},{ ia }' ^ e am plitude is then given 
explicitly by 

A({ Jab },{n)--=A^ {3ab} , {l « } )= £ lSjffc,^);^/)) f}: it (122) 

ja ja \ \ / / 
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as can be easily seen by using (|112[) and decomposing (|119|) in the SO (4) spin- network basis. This 
amplitude extends by linearity to more general states 'J = f(ip)- 

When we have a number of transitions between different 4-simplices, we have to sum over boundary 
states around each, projecting each onto this state (W\. This gives transition amplitudes for boundary 
SO(3) spin networks with the amplitude generated by the partition function 

Z =E Il( dim ^) 2 U A VfM ( 123 ) 

3f tie f v 

Observe that the quantum dynamics defined above does not change if we add the topological 
term ([75)) incorporating 7. For, by changing the integration variable to 11/ := Bf + --kBf in the 
path-integral, the above derivation goes through in exactly the same manner, and the final vertex 
amplitude differs at most by a constant, as determined by the Jacobian of the transformation from B 
to II (see appendix IC]l . 

The problem is whether this quantum dynamics defines a non-trivial theory with general relativity 
as its low energy limit. A way of systematically exploring the low-energy behavior of a background 
independent quantum theory has recently been developed (after a long search; see for instance |41p 
in [T5]. These techniques should shed some light on this problem. 

5 Conclusions 

We close with three observations. 

(i) There is a relation of the SO (4) states determined in this model to the projected spin network 
states studied by Livine in [42]. (A similar approach is developed by Alexandrov in 44 .) The 
constrained SO (4) states that form the physical Hilbert space of the model presented here can be 
constructed from (the euclidean analog of) these projected spin networks. The Euclidean analog of 
the projected spin networks defined in [42] are wavefunctions W[t/j,^ n ] depending on an SO (4) group 
element for each link, and a vector \n G SO(4)/SO(3) at each node. The wavefunctions are labelled 
by an SO(4) representation (j, + for each link, an SU(2) representation jni for each node and link 
based at that node, and an SU(2) intertwiner at each node. The SO(4) spin-networks of the present 
paper can be obtained from these projected spin networks by (i) setting j ; + = j7 = ji, (ii) setting 
jni = j; + +JT = ^ii: and (iii) averaging over Xn at each node (concretely this averaging can be done 
by acting on each Xn with an SO (4) element U n , and then averaging over the U„'s independently 
using the Haar measure). Each of these three steps corresponds directly to solving (i) the diagonal 
simplicity constraints, (ii) the off-diagonal simplicity constraints, and (iii) the Gauss constraint. 

(ii) Livine and Speziale have found an independent derivation of the vertex proposed here [45] , 
based on the use of the coherent states they have introduced in [TB] . 

(iii) The discreteness introduced by the Regge triangulation should not be confused with the 
quantum-mechanical discreteness. The last is realized by the fact that variables that give a physical 
size to a Regge cell turn out to have discrete spectrum. However, the two "discretenesses" end up to 
be related, because, due to the quantum discreteness, in the quantum theory the continuum limit of 
the Regge triangulation turns out to be substantially different that the continuum limit of, say, lattice 
QCD or classical Regge calculus. Intuitively, one can say that no further triangulation refinement 
can capture degrees of freedom below the Planck scale, since these do not exist in the theory. The 
ultraviolet continuum limit is trivial in the quantum theory (for a discussion of this point, see [3]). 
Thus, in the quantum theory the Regge triangulation receives, so to say a posteriori, the physical 
interpretation of a description of the fundamental discreteness of spacetime (the one, say, that makes 
gravitational entropy finite [46]). Numerous approaches to quantum gravity start from the assumption 
of a fundamental spacetime discreteness ([47]); others use spacetime discreteness as a regularization 
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to be removed ([48 ); the conventional formulation of loop quantum gravity derives the discretization 
of spacetime from the quantization of continuum general relativity. The derivation of loop quantum 
gravity given here is somewhat intermediate: it starts from a lattice discretization, and later finds out 
that in the quantum theory, say, the area of the triangles of the triangulation is discrete and therefore 
the regularization receives a physical meaning. See also [27] , 

In summary. We have considered the Regge discretization of general relativity. We have described 
it in terms of the Plcbanski formulation of GR. We have quantized the theory, on the basis of a 
flipped Poisson structure and imposing the off-diagonal simplicity constraints weakly. This new way 
of imposing the simplicity constraints is weaker. This weakening of the constraints is motivated by 
the observation that they do not form a closed algebra, as well as by the realization that a richer 
boundary space is needed for a correct classical limit [13] . 

The theory we have obtained is characterized by the fact that its boundary state space exactly 
matches that of (SO(3)) loop quantum gravity. This can be seen as an independent derivation of 
the loop quantum gravity kinematics, and, in particular, of the fact that geometry is quantized. A 
vertex amplitude has then been derived from the discrete action, leading to a spin-foam model giving 
transition amplitudes for loop quantum gravity states. 

We expect that the model considered here will admit a group field theory formulation [30] and 
that its vertex can be used to generate the dynamics of loop quantum gravity also in the absence of 
a fixed triangulation 3 . Whether this model is non-trivial and/or reproduces general relativity in an 
appropriate limit we do not know. Study of its semi-classical limit and n-point functions [15] should 
shed light on the discussion. 
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A 4-simplex geometry as gravitational field 

In the way we have defined the Regge triangulation above, the degrees of freedom of the gravitational 
field are captured by the geometry of each 4-simplex. Consider in each 4-simplex v Cartesian coordi- 
nates such that one vertex of the 4-simplex, say the vertex 5, has coordinates (v$) a = 0, and the other 
four vertices have coordinates (vb) a = 8?. In these coordinates, the metric will not necessarily be 6 a b, 
but rather take a form g a b{v)- The full information about the geometry of the four simplex v is then 
coded into the ten quantities g a b(v). Notice indeed that the shape of a 4-simplex in i? 4 is determined 
by 10 quantities (for instance the length of its ten sides), and therefore the ten components of g a b(v) 
are the correct number for capturing its degrees of freedom (possibly up to discrete degeneracies). 

In particular, g aa (we will often drop the argument (v) in the rest of this appendix, since we 
deal here with a single 4-simplcx) are (the square of) the lengths of the four edges adjacent to the 
5th vertex, and g a b are essentially the angles between these edges. (Also: g aa is the volume of the 
tetrahedron a, opposite to the vertex a, and g ab is the variable giving the angle between the normals 
to the tetrahedra a and b.) 

Thus, the discretized variables can be taken to be the ten quantities g a b(v). This is of course no 
surprise at all, since this is precisely a direct discretization of the variable g a b{x) used by Einstein to 
describe the gravitational field, here reinterpreted simply as a way to represent the geometry of each 
elementary 4-simplex. 

Another way of determining the geometry of a four simplex is to give its ten areas. Let A = 
1, 2, 3, 4, 5 label the five tetrahedra bounding the 4-simplex. The tetrahedron 1 is the one with vertices 
2,3,4,5, and so on cyclically. Let Jab denote the triangle bounding the two tetrahedra A and B. The 
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triangle /12 has vertices 3, 4, 5, and so on cyclically. Let Aab = Af AB be the area of the triangle Jab- 
Consider in particular the six areas A a b = Af ab , of the six faces adjacent to the vertex 5. A short 
computation, for instance using (|53[) . shows that 

A lb = 9aa9bb - gab- (124) 

Define the angle variables J aa bc = Jf ab f ac , related to the angle between the triangles f a b and f ac . (We 
can also write A\ h = J aa bb-) Then again: 

Jaabc = 9aa9bc ~ 9ab9ac- (125) 

The closure relation of the bivectors of a tetrahedron implies 

A\ 5 = Aj 2 + Aj 3 + A\ A + 2(J 1123 + Ji 134 + ^1142) (126) 

and so on cyclically. The six equations (|124[) and the the four equations (|126[) express the ten areas as 
functions of the ten components of the metric. Inverting these equations gives the metric as a function 
of the areas: g a b{AAB), and then, via (|125p . the angles as functions of the areas: Jaabc (A ab)- 
Therefore the full geometry of one 4-simplex is determined by the ten areas Aab ■ Computing the two 
functions g a b(AAB), and Jaabc (Aab) explicitly would be of great utility for quantum gravity. 

B Existence of tetrad 

In this appendix we show that the constraints (|46l l47l [48]) . at the classical discretized level, are 
sufficient to imply the existence of a tetrad associated with each tetrahedron t. The main point is 
that the 'dynamical constraint' ([49lor l50|) is not needed for the existence of e^(i). Rather, as explained 
in the main text, the role of (|49[) or (fST))) . or its reformulation as (f3"Tj) , is to ensure that the geometries 
of the tetrahedra in a single 4-simplex fit together to form a single 4-geometry. 

The triad portion of the tetrad e^(i) associated with the 3-plane of the tetrahedron t is uniquely 
determined (up to an overall sign) by the B(t)'s of the faces of the tetrahedron. The last component 
of the tetrad associated with the normal to t, however, is not determined by the B(t)'s associated 
with t. This final component is fixed only upon comparing with the tetrads at the other tetrahedra 
in the 4-simplex via the parallel transport maps Utt'(v). 

Let us make all of this precise. Let £ denote the portion of space-time associated with one of the 
4-simplices v. £ is a 4-manifold. Equip the tangent bundle of £ with a background flat connection 
(so that we have a natural notion of parallel translation between any two points, and a natural notion 
of straightness) In v, consider a tetrahedron t with faces labeled 1,2,3, and 4. For each face, we let 
B? v denote the associated bivector living in the tangent space of E. These bivectors satisfy 

4 

5>r=o. (127) 

i=l 

What we wish to show in this appendix is that the discretized constraints (|46l [47l [48|) are sufficient 
to imply that either 

Bl J = B^ei (128) 



Within a single 4-simplex, or any non-closed chain of 4-simplices, the choice of such a connection is a pure gauge 
choice. 
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or 

1 



2 

for some tetrad el(t). We begin with a 

Lemma 1. //ef, is such that B\ J = B^e^e^ for i E {1, 2, 3}, and 



*B! J := -£ /J KLBf L = B^ey v (129) 



E B " = °> ( 13 °) 



i=i 



then B\ J = B^ v e^e'l as well. An analogous statement holds with B\ J replaced by *B\ J . 
Proof. 

Using (fT30|) and then (fT27)) . 

B " = " E B " = - E S r<^ = ej#. (131) 
The proof for the *B\ J is similar. □ 

Thus, if we find a tetrad compatible with only B[ J , i?| J , B| J (in the sense of (|128[) or (|129|l ). then 
by imposing equation (|130p (the closure constraint), B\ J will then automatically be compatible with 
the same tetrad, in the same sense. In other words, by using the closure constraint, we can ignore 
B 1 /. 



Let p denote the point in the tetrahedron shared in common by faces 1,2, and 3. Let v^v^y^ 
denote the three edges emanating from this point, numbered such that face 1 is between v% and v%, 
etc, cyclically. Then we can choose our sign conventions such that 

B{ v = 2v^Vz\ and cyclically. (132) 

Specifying a tetrad is then equivalent to specifying the three internal vectors 

el :=<e£, zG {1,2,3}, (133) 

and the external unit normal to the tetrahedron. In terms of the tetrad, 

:= e^m (134) 

where nj is the unit normal to spanji^e^, v% e^, v% el}. In terms of these variables, the condition 
b ij = B^e^ei for i e {1, 2, 3} becomes 

B[ J = 2e^e^, and cyclically. (135) 
and likewise the condition *B\ J = B^ele^ becomes 

*B[ J = 2e [ ^el\ and cyclically. (136) 
We are now in a position to state the principal proposition and sketch its proof. 
Proposition 2. If B{ J , B^' 1 , B\ J are linearly independent and are such that 

*B i -B j =0, {1,2,3} (137) 

then exactly one of the two following cases hold. 
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1. There exists e\,e\,e\ such that 

B[ J = 2e^2e'^ and cyclically, (138) 

or 

2. there exists e\,e\,e\ such that 

* B[ J = 2e 2 / eg' and cyclically. (139) 

In either case, (e{, e^ef) is unique up to (e[,el,el) 1— ► (— e{, — e^, — e^). 
Proof. 

Equation (|137p . for the case i = j, using proposition (3.5.35) of [40], tells us that there exists 
<4 , PhlLH such that 

Bl J = 2aW (140) 
= 2 7 ]^/l (141) 

(i = 1, 2, 3). From the linear independence of the S's, we know in particular that none of them are 
zero. Therefore, if we define for each i 

V, := spanjc^ft} (142) 
Ui := span{7i,^}, (143) 

then Vi and Ui are each two dimensional. Vi is just the two dimensional subspace uniquely determined 
by the bivector B\ J , whereas each Ui is just the subspace uniquely determined by the bivector *B\ J . 

Let us now look at the rest of the equations (|137p . For each i ^ j, they tell us that {a i; ft, ctj, ft} 
and {7j,<5j,7j,<5j} are each linearly dependent. Vi and Vj therefore non-trivially intersect, as do U t 
and Uj, so that 

dim(Vin^) > (144) 
dim^-nUj) > 0. (145) 

But the linear independence of the B's tells us that none of the V's can be the same, and none of the 
E/'s can be the same. Thus, for i ^ j, 

dim(V-n^) < 2 (146) 

dim(Ui n Uj) < 2, (147) 

whence 

dim(Fin^) = 1 (148) 

dim(J7,- n U» = 1. (149) 

Let f( be any non-zero vector in V% n V3, etc cyclically. Likewise let f[ be any non-zero vector in 
V2 n V3, etc cyclically. One can prove that exactly one of span{/i, /a, fa} and span{/i, /2, J3} is three 
dimensional, and the other is one dimensional. To prove this is non-trivial; we leave it as an exercise 
to the reader. In proving this, one uses in a key way both the assumption that the J3's are linearly 
independent and the fact that for each i, Vi and Ui are orthogonal complements (which one can also 
show). 
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In the case where {/i,/2,/3} is linearly independent, by setting = Xifi, one can solve for A; 
such that (|138j) holds. These A, are unique up to (Ai, A2, A3) 1— > (— Ai, — A2, —A3), so that a is unique 
up to (ei,e2,e3) 1— > (— ex, — &i, — ea). Furthermore, one can show the fact dim(span{/i, /2, /3D = 1 
implies that there exists no e, such that (|139|) holds. 

In the case where dim(span{/i, fa, fs}) — 1 and /2, /s} is linearly independent, the situations 
are obviously reversed. There exist no such that (|138[) holds, whereas there exist e^'s unique up to 
(ei, e-2, ez) 1— > (— ei, — e2, —63), such that (|139p holds. 

Thus we have the proposition. □ 



C Symplectic structure 

The action is 

s = \ E ww))+5EW/ fl /)- (150) 

feintA /6<9A 

In the case when one fixes the -B variables on the boundary, one sets Rf = identity, as then the 
boundary term reduces to the usual classical boundary term appropriate when fixing the B variables 
on the boundary [36]. When one fixes the IPs (i.e., the connection) on the boundary, analogy with 
the classical theory suggests that there should be no boundary term. However, in a discrete theory 
this is not literally possible, as then no boundary variables would appear in the action at all. For 
the case of fixing the IPs on the boundary, we therefore suggest the following prescription. We want 
the terms in the boundary sum to essentially be the same as the terms in the interior sum: this 
will be the case if Rf is the holonomy around the rest of the link of / outside of A. Of course, the 
problem is that this part of the holonomy around / is not determined by any dynamical variables. 
We propose the following prescription. Whenever varying the action, keep the Rf's fixed; however, 
whenever evaluating anything on extrema of the action, set Rf equal to what we know it 'should' 
be: namely, the holonomy along the part of the link of / outside of T as determined by the BF 
equations of motion. The BF equations of motion dictate that all holonomies around closed loops 
are trivial. Thus, our prescription is, whenever evaluating an expression onshcll, set Rf = UJ 1 . A 
remark is in order as to why the BF equations of motion are used. The BF equations of motion are the 
equations of motion coming from the unconstrained action (|150p . We use the unconstrained action in 
deriving the symplectic structure because our general philosophy is that we are quantizing GR as a 
BF theory constrained at the quantum level. That is, we first compute the symplectic structure and 
quantize everything as though it were pure BF theory. Only then, after the 'kinematical' quantization 
is complete, do we impose the simplicity constraints and obtain GR. 

Let us determine the symplectic structure from the action, when we fix the LPs on the boundary. 
We will use a method essentially coming from that described in [49], and briefly mentioned in [2]; 
the method is the following. We first apply an arbitrary variation S to the action without fixing any 
variables on the boundary; then we restrict to solutions of the equations of motion. The result will be 
a boundary term 0(<5) depending linearly on the variation. This gives us a linear map from variations 
(at the space of solutions) into real numbers. This linear map is then extended off-shell basically 
by keeping the same expression in terms of the basic variables F^l We thereby obtain a linear map 
O from arbitrary variations of boundary data into real numbers — that is, a one- form on the space 



If there are constraints, this expression has ambiguities. One uses the "obvious", simplest expression in this case. 



30 



of boundary data. This is the canonical one form on the boundary phase space. The symplectic 
structure is then the exterior derivative Q — d<3. 

Let us apply this procedure. Varying the action (|150|) gives 

*S = jE Tr((5B f )U f (t)) + l ]T Tr (B f (t)(5U f (t))) 

feintA f£intA 

+ \ E Ti (( SB f) u f R f) + \ E Tr ( B f( 6u f) R f)- 

fedA fedA 

Next, restrict to the case of B's and £/'s that satisfy the equations of motion. Then, when SU / on 
the boundary is zero for / G <9A, the entire expression above must vanish. But only the fourth term 
depends on 5Uf for / S <9A. Therefore, the other terms must vanish (on solutions) no matter what 
is SU on the boundary. Therefore, on solutions, at most the fourth term above can survive, and we 
have 

Ss =lj2 ^{B f {8U f )R f ). (152) 
fedA 

Applying the equations of motion to Rf, we finally set Rf — UJ 1 1 giving us 

m) = \Y, ^(Bf(SUf)Uy 1 ). (153) 
fedA 

Extending O to arbitrary variations of the boundary data, keeping the same expression in terms of 
the basic variables, we obtain the canonical one- form 

effl^Trft^ 1 ). (154) 
fedA 

We next claim that [(SU^UJ 1 ] 1 j is an element of so(4). To show this, let Uf(A) denote a one- 
parameter path in SO(4) such that S is the tangent to the path at A = 0. Then 

(SUf)Uy 1 = Af/ /(s )| s=o{//(0 )-i = iL[t/ /(s){ / /( o)-i] s=0 (155) 

The argument in brackets is the identity at s = 0, so that one has the action of a tangent vector at 
the identity on the matrix elements of the group. Thus one has an element of the Lie algebra. 

Introduce a basis {£a}a=i....6 of so(4), orthonormal with respect to the trace — i.e. such that 
Tr(^ Q ^) = —£,a J £,pu = —5 a p. Denote the components of a Lie algebra element with respect to this 
basis just by adding a a, /?, 7 . . . superscript or subscript. (|154|) can then be written 

®(S) = ~l ^2(B f ) a [(6U f )U^] a . (156) 
fedA 

Define 

V f (6) a :=[(6U f )Uy 1 ] a (157) 

so that fij is a one- form on the copy of SO (4) parametrized by Uf. We can then write the canonical 
one form as 

9 = ~\ E (S/)«M?- (158) 
fedA 
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For each face /, and each basis element S so (4), let S^ a denote the corresponding right invariant 
vector field on the associated copy of SO (4). We then have 

N{6f, a f=€ (159) 

Thus the fij form the basis of right invariant one-forms dual to the basis S^ a of right invariant vector 
fields. These right invariant one forms have been well studied, for example, by Cartan [50]. They 
satisfy 

d/i? = ~\r P y } a m } (i6o) 

where /"/3 7 are the structure constants of the Lie algebra. Using this identity to compute the sym- 
plectic structure from 0, we obtain 



n ■■= de = ~ ( d (%)«) A M/ + \ E f a fh( B f)°i4 am}- (161) 

fedA fedA 

From this, one can read off the Poisson brackets. As there are no &B A AB terms, the 17/ 's commute. 
The first term above gives the Poisson bracket between the B's and the [7's, and the second gives the 
Poisson brackets among the -B's. Explicitly, 

{Uf,U f ,} = (162) 
{(Bf) a ,Uf,} = 25 }J ,i a U s , (163) 
{{Bf) a , (Bf)fi} = 25 fJ ,r a p(Bf) 1 . (164) 

The second two equations can be summarized by stating that the action generated by (Bf) a on Uf 
is that of the right invariant vector field field determined by 2£ Q e so (4). The above Poisson brackets 
can also be more explicitly written in terms of matrix elements 

mYj,(Uf>) K L} = (165) 
{B 1 /, {U r ) K L } = 2S fJ ,5 K V(U f )V L = 6 fjf , (S KI (U f ) J L - S KJ (U f y L ) (166) 
{{By,Bf, L } = A5 }J ,8 M ^8^ K 5 L ^ N Bf N . (167) 

If we define (t ij ) m N := 2S M ^S^ N , this becomes l|g5jl. 

Thus we obtain the symplectic structure from the action when we fix the C/'s on the boundary. The 
case of fixing the B's on the boundary is not as well-understood. Application of the above prescription 
to this case seems to lead to a non-SO(4)-gauge-invariant symplectic structure in which all the B's 
commute. For the present, we simply do not address this problem, and take the symplectic structure 
to be the one determined with ?7's fixed. 

So far we have seen that the action we have considered gives the unflipped Poisson structure. 

Recall, however, that in the LQG approach the action that is quantized is the GR one with a 
topological term that doesn't change the equations of motion: 

Slqg = J *(eAe)AF+iy (eAe)Af (168) 

where 7 is the so called Immirzi parameter. Recall that the introduction of this topological term is 
required in order to have a theory of connections on the boundary: without it, as shown by Ashtekar, 
the connection variable does not survive the Legendre transform [34]. Let us therefore consider the 
discretization of the modified BF action: 



S= I [B + -*B ) A F + 4>uklB 1j A B kl (169) 
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and follow the lines presented in this paper. Define then 11/ :— Bf + A * Bf, the symplectic structure 
is then such that, for each face, the Uy are identified with the generators of SO(4). The constraints 
Cff imply, in terms of these new variables: 



Note that for 7 <C 1 and 7> 1 one recovers the same constraint *TLfH.f ~ being used before, whence, 
in these two limits, the diagonal simplicity constraints again imply simplicity of the representations of 
the edges (j + = (Thanks to Alejandro Perez for pointing this out.) Furthermore, as mentioned 
in H4.31 neither does the dynamics change with the introduction of the 7-term in the action. What 
is different, then, with the introduction of the 7 term, when 7 < 1 or 7 > 1? In the 7 <C 1 case, 
the symplectic structure in terms of B is different — it is the flipped symplectic structure. Is there 
an independent reason to prefer 7 <§; 1 over 7> 1? Yes: the latter is unphysical, because it leads to 
a macroscopic area spectrum in the boundary theory. Thus, when we look at things more carefully, 
we do in fact see that the flipped symplectic structure is a consistent choice. And with this flipped 
symplectic structure, as noted in the main text, the off-diagonal simplicity constraints, in the form 
(|105p . lead precisely to the model presented in this paper. 

As a side remark, consider the case 7> 1. In this case, one obtains the more standard symplectic 
structure ([86]) . As noted in §4.2[ this symplectic structure, along with the constraint (|105p . imposed 
as discussed above, leads to the Barrett-Crane model. Thus the 7 ^> 1 case, which we argue against, 
can be viewed as yielding the Barrett-Crane model. 
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